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Understanding a behavior of galaxy biasing is crucial for future galaxy redshift surveys. One aim 
is to measure the baryon acoustic oscillations (BAOs) within the precision of a few percent level. 
Using 30 large cosmological iV-body simulations for a standard ACDM cosmology, we study the halo 
biasing over a wide redshift range. We compare the simulation results with theoretical predictions 
proposed by Matsubara which naturally incorporate the halo bias and redshift-space distortions into 
their formalism of perturbation theory with a resummation technique via the Lagrangian picture. 
The power spectrum and correlation function of halos obtained from Lagrangian resummation theory 
(LRT) well agree with TV-body simulation results on scales of BAOs. Especially nonlinear effects on 
the baryon acoustic peak of the halo correlation function are accurately explained both in real and 
redshift space. We find that nonlinearity and scale dependence of bias are fairly well reproduced by 
1-loop LRT up to k — 0.35/iMpc -1 (z = 2 and 3) within a few percent level in real space and up to 
k = O.lfrMpc -1 (z = 2) and 0.15/iMpc" 1 (z =3) in redshift space. Thus, the LRT is very powerful 
for accurately extracting cosmological information in upcoming high redshift BAO surveys. 
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PACS numbers: 98.80.Es 



I. INTRODUCTION 



The redshift survey of galaxies is one of the most im- 
portant methods to probe the nature of a mysterious en- 
ergy component called dark energy, which is supposed 
to explain the late-time cosmic acceleration discovered 
by the observation of distant supernovae [l|, • Baryon 
acoustic oscillations (BAOs) imprinted on the large scale 
structure can be used as a standard ruler to measure 
the cosmic expansion history of the Universe [e.g., 0-0] ■ 
The characteristic scale of BAOs, which is determined 
by the sound horizon scale of baryon-photon plasma at 
the recombination epoch, is though to be a robust mea- 
sure [e.g., HHnJ. Using the BAOs as a standard ruler, 
large galaxy surveys are expected to provide a tight con- 
straint on the nature of dark energy [e.g.. [TT| - [T6j |. 



That is why most of the planned galaxy redshift sur- 
veys aim at measuring BAOs within the precision of a 
few percent level. Some examples are BigBOSS (l7j . 
Euclid [l8[, the Hobby-Eberly Telescope Dark Energy 
Experiment (HETDEX) 1 , the Large Sky Area Multi- 
Object Fiber Spectroscopic Telescope (LAMOST) 2 , the 
Square Kilometre Array (SKA) 3 , the Subaru Measure- 
ment of Images and Redshifts (SuMIRe) 4 and the Wide- 
Field Infrared Survey Telescope (WFIRST) 5 . Recent ob- 
servations of BAOs in modern galaxy surveys work well 
and constrain the cosmological parameters with approx- 
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1 http: / /hetdex.org/ 

2 http:/ /www. lamost.org/ website/en 

3 http://www.skatelescope.org/ 

4 http://sumire.ipmu.jp/en/ 

5 http://wfirst.gsfc.nasa.gov/ 



imately 10% level [see, e.g.,ll9ri2 

However, in order to pursue an order-of- magnitude 
improvement, a theoretically precise description of the 
BAOs is a crucial issue. It needs to be investigated taking 
account of the various systematic effects. The two main 
systematic effects are galaxy biasing and redshift-space 
distortions. The spatial pattern of galaxy distribution is 
not the same as that of dark matter and the galaxies are 
biased tracers of mass 0, [13] . The redshift of a galaxy 
does not purely reflect the Hubble flow, and Doppler shift 
by a peculiar velocity is inevitably added. Therefore the 
spatial pattern of clustering of galaxies is distorted in 
redshift space [28f . With time, these effects are increas- 
ingly difficult to treat due to nonlinear evolution of struc- 
ture growth and degrade the contrast of the BAOs in the 
matter power spectrum and/or correlation function, de- 
creasing the signal-to-noise ratio of the standard ruler 
method. The nonlinear effects on BAOs play an impor- 
tant role in taking into account percent-level cosmology 
although baryonic features appear on large scales such 
as 100ft. - : Mpc. The amplitude of BAOs is increasingly 
damped from small scales to large scales, with decreasing 
redshift [e.g.. [29l - [37j . Redshift distortion effects enhance 
a nonlinear damping of BAOs along the line-of-sight di- 
rection [e.g.. |29W37| . Biasing is also affected by nonlin- 
ear effects and is scale dependent on BAOs scales. Such 
nonlinear effects im pos e a serious problem in analyzing 
galaxy surveys [e.g., TH UlEj. 

Since galaxies are expected to form within dark matter 
halos, understanding and modeling the clustering prop- 
erties of the halos play an important role in model galaxy 
biasing. In the usual halo model approach, the halo clus- 
tering is modeled by linear dynamics and linear bias fac- 
tors [4^445j . Recently, Smith et al. [HHil and Elia et aL 
[47j developed a nonlinear perturbation theory, incorpo- 
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rating the effects of halo bias. In their approach, there 
are some problems. Smith et al. [39|, H(| treat the halo 
bias as a local bias in Eulerian space, although the halo 
bias is intrinsically nonlocal in Eulerian space and does 
not fit well into the local Eulerian biasing scheme [48| . 
Meanwhile, the analytic approach developed by Elia et al. 
[47| cannot deal with the redshift-space distortions. 

In this paper, we use nonlinear perturbation theory de- 
veloped by Matsubara 49] which naturally incorporate 
the halo bias and redshift-space distortions in their for- 
malism of perturbation theory with a resummation tech- 
nique through the Lagrangian picture. A significant ad- 
vantage of the Lagrangian resummation theory (LRT) 
is that it is simpler and easier to calculate the power 
spectrum than other resummation methods even in the 
presence of halo bias and redshift-space distortions. The 
computational cost is similar to that of standard pertur- 
bation theory (SPT) [e.g.,[5(|. In this work, we exam- 
ine how well it reproduces the simulation results in both 
real and redshift space. We focus not only on the halo 
power spectrum but also on the halo correlation function, 
because cosmological information that can be extracted 
from them is not exactly equivalent to each other because 
of different error properties. In this paper, we use 1-loop 
LRT for simplicity. The 2-loop corrections should be use- 
ful to extend the valid range, once we confirm that the 
1-loop results agree with simulations [5l| . 

The structure of this paper is as follows. In Section. [TTl 
we briefly review analytical models in order to obtain 
matter power spectrum in real and redshift space. We 
also review the halo power spectrum based on 1-loop 
LRT. In Section. Mil we describe the details of TV-body 
simulation and a method to calculate the power spectrum 
and two-point correlation function from TV-body simula- 
tions. In Section. IIVI we show the main results. SectionlVl 
is devoted to our conclusions. 



II. ANALYTICAL MODELS 

There are several analytical models to account for evo- 
lutions of the matter power spectra in real and redshift 
space. In this paper, we mainly use three analytical mod- 
els to compare the TV-body simulation results; Linear per- 
turbation theory (LIN) [e.g., [jj, 1-loop SPT [e.g.. [50j. 
and 1-loop LRT [e.g., H3Ej. 



A. Matter power spectra in real space 

Assuming that the amplitude of fluctuation is small, 
we can derive SPT as expansion of the fluid equations. 
Schematically, the expansion is written as 

P SPT (k, z) = P+(z) 2 P L , (fc) + D+(z)*P^(fc) + • • • , 

(1) 

where D + (z) is the linear growth factor normalized to 
unity at present and -Pl.o(^) an d Pspto'W are the 



linear power spectrum and the 1-loop contribution to 
power spectrum, respectively, at present. For the sake 
of convenience, by using Pj J (k,z) = D + {z) 2 Pi J fi{k) and 
Pgp™ p (fc, z) = D + (z) 4 Pgp^°Q (k), we rewrite Eq. ([1]) as 



PsPT(fc)=P L (fc)+P s 1 p 1 ° OP (fc) 



(2) 



Here we drop the z dependence. Throughout this paper, 
we do not write the z dependence unless otherwise stated. 
The explicit expressions of the 1-loop contribution can be 
found in the literature [e.g.. [62i - [66| . 

On the other hand, 1-loop LRT is written as [6l| 

P L R T (fc) = exp (-fcV 2 ) [P SPT (fc) + k 2 alP L (k)] , (3) 

in terms of SPT. Here a 2 is the one-dimensional linear 
velocity dispersion given by 



d 3 p Pl(p) 
(2tt) 3 p 2 ' 



(4) 



If the exponential prefactor is expanded, the 1-loop SPT 
is exactly recovered when we consider the 1-loop level. 

Finally we comment on other perturbation theories. 
Several other approaches have been proposed beyond 
SPT, such as the renormalized perturbation theory 
(RPT) [13-[69|L the large- N expansion [z3|, the Time- 
RG method [7ljj , the renormalization group approach [72[ 
and the closure theory [73[ . In these newly developed ap- 
proaches, the standard perturbative expansion is reorga- 
nized and partially resummed in various ways. Different 
levels of approximations and ansatz are used in those ap- 
proaches. 



B. Matter power spectra in redshift space 

As in real space, we use SPT and LRT in redshift space. 
Schematically, SPT in redshift space is written as 

P s S p T (fc,^ = (l + ^ 2 ) 2 PL(fc) + Ps S p 1 T OOP (fc>M) + --- , (5) 

where / = d\n.D/d\na is the logarithmic derivative of 
the linear growth rate and \i is the cosine of the angle 
between the line of sight and the wave vector k. Here 
PgPT loop (fc) is the 1-loop contribution to the redshift- 
space power spectrum in SPT [6l[ • The first term on the 
right-hand side is the linear-order result of the redshift- 
space power spectrum and the factor (1 + ffJ 2 ) 2 indicates 
the enhancement of the power spectrum which is the so- 
called Kaiser effect [28| . The Kaiser effect is represented 
as the coherent distortion by a peculiar velocity along the 
line-of-sight direction. 

On the other hand, 1-loop LRT in redshift space is 
written as [6l[ 

Ptarfcn) = exp [-k 2 a 2 {l + /(/ + 2) f i 2 }] 



(i + V)^p L (fc) + p s ™M 



+ (1 + /^ 2 ) 2 [1 + /(/ + 2)/i 2 ]A; 2 a 2 P L (fc) 



(6) 
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As in the case in real space, if the exponential prefactor 
is expanded, Eq. © reduces to the 1-loop SPT result. 
The exponential prefactor shows the nonlinear damping 
effect by the random motion of peculiar velocities, i.e., 
the Fingers-of-God effect [zllzi- 

We also use a model of redshift-space power spectrum, 
proposed by Scoccimarro [76| . This model gives 



*sco(M) 



exp (-/VfcX) 

[P u (k) + 2ffi 2 P 5e (k) + f^Pge] , (7) 



where P$s(k), Pgg(k) and Pse(k) are auto power spectra 
of density and velocity, and their cross power spectrum, 
respectively. This model accounts for the nonlinear ef- 
fects although it is still phenomenological. 

In this paper, we consider only the monopole power 
spectrum which is defined as 



P s (k) 



1 



d^P s (k,n). 



(8) 



Other multipole contributions to the power spectrum 
and correlation function in redshift space, such as the 
quadrupole and hexadecapole are examined in the liter- 
ature [e.g., [77" 



C. Halo power spectra in real and redshift space 

Based on 1-loop LRT, the power spectrum of the biased 
object in redshift space is given by 49] 

^hh,LR T (fc,M) = ex P [- fc ^v[l + /(/ + 2)M 2 ]] [(1 + (F') 



f^) 2 P L (k) + J2^ n f m E nm (k)}, (9) 



where (F') denotes a Lagrangian linear bias factor and 
the explicit expression of E nm (k) is shown in Matsubara 
[49l | . E nm (k) includes the higher-order Lagrangian bias 
factor such as (F"). Up to 1-loop order, we need only 
two bias factors (F') and (F"). When we consider a mass 
range [Mi, M2] of halos, the bias function for the halo 
bias is given by [iij ] 



I. - J-J J Mi v du" 



dM M 



r M 2 / \ 



d In a dM 
dM M 



(10) 



where S c is the critical overdensity at the present time. 
In an Einstein-de Sitter cosmology, S c — 1.686. For gen- 
eral cosmology, the value of S c shows weak dependence 
on cosmology [86H89| | , so we include cosmological depen- 
dence on S c . The quantity a is the root-mean-square 
linear density smoothed with a top hat filter of radius R 
and enclosing an average mass M = pqAttR 3 /3, 



a 2 (M) 



k 2 dk 



W 2 (kR)P L (k), 



(11) 



with 



W{x) 



-(sin a; — xcosx), 



(12) 



where po is the mean matter density of the Universe. 
The quantity v is defined by v = S c /a. g(a) is the scaled 
differential mass function defined as 



M dn 
po diner" 



(13) 



where n is the number density of halos with mass M. 
The quantity g(a) is frequently used in the literature and 
there have been several anal ytic al predictions [9TM93I ] and 
fitting formulae [e.g., Ill [111)1] for g(a). In section[lVC] 
we will compare the mass function obtained from our N- 
body simulations with several fitting formulae to examine 
which fitting formulae are better. 

From Eq. © with substituting / = 0, the power spec- 
trum of the biased object based on 1-loop LRT in real 
space is obtained as 



hh.LRT 



(k) = cxp [-k 2 o 2 ] [(1 + (F')) 2 P L (k) + E 00 (k)] . 

(14) 



By expanding the exponential prefactor and considering 
the linear term in I\,(k), the linear result is derived as 



Phh,LINW = (l + (F')) 2 PL(fc). 



(15) 



In linear theory, using the Eulerian linear bias factor b, 
the power spectrum of the biased object is defined by 



Phh,Lw(k) — b^PjjQz). 



(16) 



Comparing Eq. (IT5|) and Eq. ([T6]) . we can easily derive 
the relation between Eulerian and Lagrangian biases as 



6=1 + (F') 



(17) 



This result is the same as that derived from the halo 
mo del a pproach by using the spherical collapse model 
[42l I lOOt ] . i.e., the Eulerian linear bias is given by unity 
plus the Lagrangian linear bias (see, [48] for more accu- 
rate expressions on nonlinear scales). It should be noted 
that the result of Eq. (fl7|) is derived without assuming a 
spherical collapse model. 

The corresponding linear result in redshift space is ob- 
tained by a linear term of Eq. (|9|) as 



P 



hh.LIN 



(k,v) = (l + (F>)+fv 2 ) 2 P L (k). (18) 



By using Eq. (|17[) . this equation is rewritten by 

PLxw(k,v) = b 2 (l + /3p 2 ) 2 P L (k), (19) 

where f3 = f/bis the redshift-space distortion parameter. 
This equation is equivalent to the Kaiser formula [28| . 
Methods used to determine the redshift-space distortion 
parameter ft are summarized by Hamilton [771 ] . 
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TABLE I: Parameters used in low- and high-resolution iV-body simulations. 



Name 




Oa 


fib 


h 


n s 


08 




N p 


Zini 


r s 


Ntun 


L1000(low resolution) 


0.265 


0.735 


0.0448 


0.71 


0.963 


0.80 


1000/1 


_1 Mpc 


1024 3 


36 


50/i -1 kpc 


30 


L500(high resolution) 


0.265 


0.735 


0.0448 


0.71 


0.963 


0.80 


500ft- 


x Mpc 


1024 3 


42 


25ft~ 1 kpc 


5 



III. iV-BODY SIMULATIONS 
A. Simulation Parameters 

We use the numerical simulation code Gadget2 [lOlj 
in its full tree-particle mesh mode. For all of the sim- 
ulations discussed in the present paper, we adopted 
the standard ACDM model with the matter density 
Cl m — 0.265, baryon density f2& = 0.0448, dark en- 
ergy density Oa = 0.735 with equation of state pa- 
rameter w = — 1, spectral index n s — 0.963, the vari- 
ance of the density fluctuation in a sphere of radius 8 
/i _1 Mpc (78 = 0.80, and Hubble parameter h — 0.71. 
These cosmological parameters are consistent with the 
Wilkins on M icrowave Anisotropy Probe 7yr (WMAP7) 
results |l02j |. We employ N p — 1024 3 particles in boxes 
of side ibox = 1000/i _1 Mpc with softening length r s be- 
ing 50/i _1 kpc and of side Lbox = 500/i _1 Mpc with soften- 
ing length r s being 25/i _1 kpc, respectively abbreviated to 
L1000 and L500. Throughout this paper, we show results 
obtained from L1000 unless otherwise stated. The ini- 
tial conditions are generated based on the 2nd- order La- 
grangian perturbation theory (2LPT) [l?53l[l04T| with the 
initial linear power spectrum calculated by CAMB [105J. 
We use parallelized 2LPT code which is kindly provided 
by Takahiro Nishimichi who developed it in Valageas and 
Nishimichi jl06| to run large cosmological A^-body simu- 
lations with initial conditions based on 2LPT. The initial 
redshift is set to z- m i = 36 for L1000 and Zi n i = 42 for 
L500. We perform N vun = 30 realizations for L1000 and 
N lna = 5 realizations for L500. Table Q] summarizes the 
parameters used in the simulations. 

We store outputs at z = 3, 2, 1, 0.5, 0.3, and and 
identify halos for each output using a friends-of-friends 
(FOF) group fin der w ith linking length of 0.2 times the 
mean separation jl07| . We select halos in which the num- 
ber of particles, N p , is equal to or larger than 20 which 
corresponds to the halos with masses 1.37 x 10 12 /i _1 M© 
for L1000 and 1.71 x lO u ft _1 M for L500. The average 
number, number densities, and mass of halos among re- 
alizations for redshifts at which we store outputs can be 
found in Table ED 



B. Analysis: Power Spectra and Two-Point 
Correlation Functions 

To calculate the power spectrum from A^-body simu- 
lations, we calculate the Fourier transform of the den- 



sity field, denoted as 5 n (k), where the superscript n de- 
notes the n-th realization and k shows the wave number. 
First, we assign the A^-body particles onto N^ rid = 1024 3 
grids ba sed on the cloud-in-cell (CIC) mass assignment 
scheme [108| . We then use Fast Fourier Transformation 
(FFT) to calculate the density contrast in Fourier space 
an d correct t he effect of the CIC mass assignment scheme 

as pi [not 



6 n (k) -> 



S n (k) 



sine 



2 ' 



. ( 20 ) 

where sinc(x) = si n(x)/a:. A more careful analysis was 
done by Jing [lll| . Finally, squaring the density con- 
trast in Fourier space and taking an average over Fourier 
modes and realizations, the ensemble average of binned 
power spectrum is given by 



JV rnn 



1 



E 

kf In <\k\<k™ 



»<ifci<fc? 



S n (k) 



(21) 
(22) 



where Nf is the number of Fourier modes in the i-th wave 
number bin, and fcf" n and fc™ ax are the minimum and 
maximum wave number of the i-th bin, respectively. The 
obtained power spectrum is contaminated by the shot 
noise effect due to the discreteness of density field in A^- 
body simulations. We simply assume the Poisson model 
where the shot noise is given by the inverse of number 
density of dark matter particles, 1/n = L^ oy ./N p and then 
subtract the shot noise, though this effect is very small 
and only has an impact on the result at z — 3 on scales we 
have considered. Note that the power spectra measured 
from the treatment above suffer from the effect of finite- 
mode sampling due to a finite number of real izati ons or 
a finite box size, shown by Takahashi et al. [l!2j . The 
finite-mode effect is important at k 0.1/iMpc -1 . We do 
not eliminate finite- mode effect because extension of their 
method to the halo case we are interested in is not trivial. 
Hence, to overcome this effect we run many realizations, 
i.e., A^run = 30 and use a large simulation box size, i.e., 
L box = 1000/i" :L Mpc. 

To calculate the two-point correlation function from 
Af-body simulations, we adopt the grid-based calculation 
using FFT. Taruya et al. [73] show a grid-based calcu- 
lation with FFT almost coincides with the direct pair- 
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TABLE II: Properties of halo catalogues for low- and high-resolution simulations. Nh, rn and Mu are the ensemble average 
number, number densities and mass of halos for redshifts we employed. 



L1000(low resolution) 


L500(high resolution) 






N h 


n h [h 3 Mpc- 3 ] 








N h 


n h [h 3 Mpc- 3 ] 




Z 


= 3 


4.00x10 s 


4.00xl0" 4 


2.59xl0 12 


z 


= 3 


1.08xl0 6 


8.68xl0 -3 


4.58xlO n 


Z 


= 2 


1.21xlO B 


1.21xl0" 3 


3.30xl0 12 


z 


= 2 


1.86xl0 6 


1.48xl0" 2 


6.12xlO n 


Z 


= 1 


2.38xl0 6 


2.38xl0 -3 


4.75xl0 12 


z 


= 1 


2.42xl0 6 


1.94xl0 -2 


9.07x10" 


Z 


= 0.5 


2.82xl0 6 


2.82xl0" 3 


5.99xl0 12 


z 


= 0.5 


2.52xl0 6 


2.01xl0" 2 


1.15xl0 12 


Z 


= 0.3 


2.93xlO B 


2.93xl0" 3 


6.63xl0 12 


z 


= 0.3 


2.52xl0 6 


2.01xl0" 2 


1.27xl0 12 


Z 


= 


3.05xl0 6 


3.05xl0" 3 


7.73xl0 12 


z 


= 


2.49 xlO 6 


1.99xl0 -2 


1.47xl0 12 



count method (see, Appendix C in their paper). Since 
the grid-based calculation with FFT is computationally 
less expensive than the direct pair counting, we therefore 
adopt the grid-based calculation. It should be noted that 
the grid-based calculation is limited to scales larger than 
the grid size r > ^box/^Vgrid- In this method, we first 
compute the square of the density field in Fourier space 
on each grid. Then taking the inverse Fourier transfor- 
mation and an average over distances and realizations, 
we can obtain the two-point correlation functions. This 
is expressed as [79| 



Un) = 



l 



N„ m NJ 



E 

n=l 



E 



FFT" 



5 n (k) 



(23) 

where FFT 1 denotes the inverse FFT of the square of 
the density field in Fourier space. Here, N[ is the number 
of modes in the i-th distance bin, and rf 1 m and rf ax are 
the minimum and maximum distances of the i-th bin, 
respectively. We chose to be the center of the i-th bin, 
i.e., n = (r™ in + rf iax )/2. 

For the estimation of the halo power spectrum, we 
adopt the same method used for the power spectrum es- 
timation to the halo cases, but apply the FOF halo mass 
correction as will be discussed in Section. IIVCI To be 
more precise, we impose Eq. (|29p to halo particles when 
we assign the halo particles based on the CIC mass as- 
signment scheme in order to introduce the FOF mass cor- 
rection for A^-body simulations. Then replacing 8 n with 
6%, we do the same calculation described above [from 
Eqs. (|2"0)) - (P22"l) ] in order to obtain the halo power spec- 
trum. 6% is the density field of halos in Fourier space in 
the n-th realization. 

There is one concern about shot noise. If the dark 
matter halos are regarded as a Poisson sampling, the 
subtraction term is the same as the dark matter particle 
case but using the appropriate number density n^- How- 
ever Smith et al. [39| found that this standard correction 
method is not exactly correct for halos, particularly for 
those of large mass. This is probably because in order to 
identify halos by using the FOF algorithm, we automati- 
cally impose that distances between halos are larger than 



the sum of their radii, or they would have been linked as 
bigger halos. To correctly eliminate the shot noise, Smith 
et al. [3!| proposed a new procedure that includes this 
exclusion effect (see, Appendix A in thei r paper for more 
details and also see, Seljak et al. (lT3 |, Hamaus et al. 



114] for a new method to suppress the shot noise effect 
by weighting halos). However, we do not use this pro- 
cedure to subtract the shot noise, but use the standard 
correction method for simplicity. The result could not 
being changed on scales we are interested in, because the 
number density of halos are large enough, and there is 
little shot noise effect. 

For the estimation of the halo two-point correlation 
function, we use a grid-based calculation with FFT as is 
used in matter correlation function. 



IV. RESULTS 
A. Matter Power Spectra 

FigureQ]shows the matter power spectra obtained from 
A^-body simulations among 30 realizations for z = 3, 2, 
1 and 0.5. We compare them with several analytical pre- 
dictions. The left and right panels show the results of 
power spectrum divided by the smooth linear power spec- 
trum P nw (k) which is calculated by u sing the no-wiggle 
fitting formula of Eisenstein and Hu [ll5| in real space 
and by the no-wiggle power spectrum taking the Kaiser 
effect (1 + 2//3 + / 2 /5) into account in redshift space, 
respectively. 

All of the theoretical predictions well agree with the N- 
body simulations within a range of the error bars on large 
scales. The range of agreement in all of the theoretical 
predictions is generally wider as redshift is higher. This 
is because the amplitude of density fluctuation is smaller 
at highe r redshift. As has been stated in the literature 
[IE Il04| , the result of 1-loop SPT (dashed curves) is not 
sufficiently accurate to describe the BAOs not only in real 
space but also in redshift space. It seems that the discrep- 
ancy in redshift space between the 1-loop SPT prediction 
and simulation results is larger than that in real space. 
The amplitude of the power spectrum from 1-loop SPT 
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FIG. 1: Comparison of mass power spectra obtained from iV-body simulations to analytical predictions in real (left panel) and 
redshift (right panel) s pace for redshifts, z = 3, 2, 1, and 0.5 from top to bottom. Power spectra are normalized by the no-wiggle 
fitting formula P n -w(k) [115| ] in real space and by the no- wiggle fitting formula taking the Kaiser effect into account P^ w (fe) in 
redshift space, respectively. The solid, dotted and dashed curves show the predictions of the 1-loop Lagrangian resummation 
theory (1-loop LRT), the linear theory (LIN), and the 1-loop standard perturbation theory (1-loop SPT), respectively. The 
dot-dashed curves in the right panel show the predictions of the empirical model proposed by Scoccimarro [7fj]. The arrows 
indicate the valid ranges for several redshifts, where the results of 1-loop LRT are expected to be accurate within a few percent. 



generally overestimates that from iV-body simulations. 
However, the range of the agreement in 1-loop SPT is 
wider in both real and redshift space than that in LIN 
(dotted curves), because 1-loop SPT includes the next- 
order contribution of nonlinear growth. Meanwhile, the 
amplitude of the power spectrum in 1-loop LRT (solid 
curves) rapidly falls off at a certain wave number and 
it deviates from the iV-body results. This is attributed 
to the exponential prefactor in Eqs. ([3]) and ©. In real 
space, the agreement of 1-loop LRT is roughly equiva- 
lent to that of 1-loop SPT. However, we can see that the 
1-loop LRT results give a better agreement with simula- 
tions in redshift space. We also plot another theoretical 
prediction proposed by Scoccimarro [7o| in redshift space 
which is denoted as the dot-dashed curves in the right 
panel of Figure [TJ The model of Scoccimarro [76j seems 
to be a good agreement at certain redshift (z — 0.5). 
However, this is just coincidence because it does not de- 
scribe the simulation results at different redshifts. 

Comparing the power spectra in real and redshift 
space, the amplitude of the power spectra on small scales 
in redshift space is suppressed by nonlinear redshift-space 
distortions. This suppression is due to the random mo- 
tion of peculiar velocities in virialized objects which is 
known as the Fingers- of-God effect 0, [Zg. 

The arrows in Figure [T] indicate the reliable ranges 
&:nl/2 for several redshifts, where the results of 1-loop 
LRT are expected to be accurate within a few percent. 
This criterion is proposed by Matsubara 61]. It is de- 



termined by the damping scale of the exponent which is 
estimated from Eq. (|4]) as 



1 



6vr 2 



d P P h (p) 



1-1/2 



(24) 



B. Correlation Functions 

Next we focus on the two-point correlation function. 
Although the correlation function and power spectrum 
are directly related by Fourier transforms and have math- 
ematically equivalent information, cosmological informa- 
tion that can be extracted from them with real data is 
not exactly equivalent because error properties are dif- 
ferent. Therefore it is important to examine not only the 
power spectrum but also the two-point correlation func- 
tion. Two-point correlation functions in real and redshift 
space around baryon acoustic peaks at z — 3, 2, 1, 0.5, 
and are shown in left and right panels of Figure [2j re- 
spectively. The theoretical two-point correlation function 
can be expressed in terms of the power spectrum as 



k 2 dk sin(fcr) 



2tt 2 



kr 



P(k). 



(25) 



The ./V-body results clearly deviate from the LIN predic- 
tion (dotted curves) as decreasing the redshift, because 
nonlinear growth of the structure becomes significant as 
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FIG. 2: Comparison of two-point correlation functions obtained from JV-body simulations to analytical predictions in real (left 
panel) and redshift (right panel) space for redshifts, z = 3, 2, 1, 0.5, and from bottom to top. The solid curves represent the 
results of the 1-loop Lagrangian resummation theory while the dotted curves show the results of the linear theory. 



decreasing the redshift. In contrast, 1-loop LRT results 
(solid curves) fairly well reproduce the TV-body results 
at all the redshift we have considered in both real and 
redshift space, although in redshift space a closer look 
around acoustic peaks reveals that the iV-body results 
deviate from 1-loop LRT results as decreasing the red- 
shift. Although the 1-loop LRT predictions considerably 
deviate from A-body simulation results on large wave 
number in the power spectrum (Figure [TJj the acous- 
tic peak structure in the correlation function comes from 
the low-fc behavior of the power spectrum, and the power 
spectrum at low-A: is accurately described by 1-loop LRT 

[sum. 

The baryon acoustic peaks in both real and redshift 
space tend to be smeared as decreasing the redshift. 
However, the effect seems stronger than those in real 
space due to the Fingers-of-God effect. In the power 
spectrum results (Figure [1]), the difference between 1- 
loop SPT and 1-loop LRT seems to be small, especially 
in real space. In correlation function results, however, 
the difference is quite big, because SPT cannot predict 
the correlation function which is ascribed to fail to con- 
vergence the integral in Eq. (|2"5j) [6lj . Therefore, we con- 
clude that the LRT is very powerful and sufficient to 
predict the BAOs in correlation functions. 



C. Halo Mass Functions 

If we want to accurately calculate the halo power spec- 
tra based on LRT, accurate halo mass functions are 
clearly needed, since the Lagrangian bias factor is ex- 



pressed in terms of the halo mass function (see, Eq. [TO]) . 
Here we use different mass function expressions given pre- 
viously to compare our A-body simulation results. A fa- 
mous numerical fit for g(a) is given by Sheth and Tormen 
[94| (hereafter ST) which is expressed as 



1 



aS*J . 



exp 



aSl 
2a 2 



(26) 



with A = 0.3222, a = 0.707 and p = 0.3. Recently, 
Crocce et al. [l?} (hereafter MICE) recalibrated the halo 
mass function using a large set of A-body simulations 
with good mass resolution and large cosmological vol- 
umes called MICE simulations. They provided a numer- 
ical fit as 

c(zy 

rr2 



5mice(ct) = A(z) 



-a(z) 



+ b(z) 



exp 



(27) 



-0.15 



with A{z) = 0.58(1 + z)~ - 13 , a(z) = 1.37(1 + z) 
b(z) = 0.3(1 + z)-°- 084 , and c(z) = 1.036(1 + z)-°- 024 . 
Bhattacharya et al. [99j (hereafter Coyote) also recently 
recalibrated the mass function using their simulations so 
called as Coyote simulations whose parameters are very 
close to MICE simulations. Therefore, both results show 
good agreement within a few percent except at very high 
masses. The functional form of the Coyote fit which is 
similar to the ST mass function is expressed as 



5 , Coyotc(c r ) = A\J - 



1 



a 

m 



exp 



2<7 2 



(28) 



8 



o 
ft 



10 11 r 



10 1 



S 
ao 
_o 

\ 
c 

T3 



10 9 



10 8 



i i 1 1 in — i i 1 1 1 i — i i 1 1 mi — i 1 1 1 in — r~rr 
1 _ # t 



-v 

no correction n«U 
(offset factor 10) "^"^ 



10 7 
10 11 




■ L1000 at z=0 

= L500 at z = 

.• L1000 at z=l 

o L500 at z=l 

— Bhattacharya et al. 

Sheth & Tormen 
I I i 



10 



12 



10 13 10 14 
M h [h"> M Q ] 



10 1 



o 



s 

M 
O 



T3 



10 1 



10 1 



I III 1 III III 1 1 1 1 1 I I I I III 1 Mill 



no correction 

(offset factor 10)?% ~> 




10 la 10 13 10 14 

M h [h-' M Q ] 



10 1 



FIG. 3: Halo mass function at redshift 0, 1 (left panel), 2 and 3 (right panel). We calculate the halo mass function from 
iV-body simulations both with the force plus N p correction and without the correction (see text for correction in the details). 
The latter was shifted upward by a factor of 10 for clarity. The solid and dashed curves denote Bhattacharya et al. [9!j and 
Sheth and Tormen [§3] results, respectively. 



with A = 0.333(1 + z)-°- n , 5 = 0.788(1 + z)" 01 , p = 
0.807, and q = 1.795. 

In Figure [3J we plot the halo mass function as a func- 
tion of halo mass at z = 0, 1, 2, and 3. We show the 
results of L1000 as filled symbols and the results of L500 
as open symbols, respectively. The simulation results are 
plotted with the force plus the N p correction proposed by 
Bhattacharya et al. [99( and without the correction. The 
latter was shifted upwards by a factor of 10 for visibility. 
The net correction for the FOF halo mass is expressed as 




M c JM h = [1.0 - 0.04(r s /650kpc)](l - n" ' 65 ), (29) 

where Mfc is a corrected mass, Mh is an uncorrected 
mass, and nn is the number of particles that construct a 
halos. This equation is slightly modified from the origi- 
nally suggested equation by Warren et al. [95| : 

M c h /Mh = 1 - n- - 60 . (30) 

The solid and dashed curves denote the results of Coy- 
ote and ST fitting formulae. Our simulation results well 
reproduce the Coyote results over a wide range of halo 
masses and redshifts, while ST slightly deviates from the 
simulation results, especially at large halo masses. The 
impact of the correction to the FOF mass can be easily 
seen in Figures 0] and 

In Figure 31 to investigate the agreement in more quan- 
titative ways, we show the ratio of the halo mass func- 
tion obtained from iV-body simulations to that from fit- 
ting formulae of Coyote and ST used in Figure El In 
addition to the above two fitting formulae, we use the 



recently proposed fitting formula by MICE. The left and 
right panels show the results of L500 and L1000 at z = 
and 1. We can see good agreement between our sim- 
ulation results and the Coyote fitting formula within a 
few percent if we apply the correction (Eq. [29]) to FOF 
halo masses. Also, the MICE fitting formula reproduces 
the simulation results applying FOF mass correction with 
the same level in Coyote at z = 0, although the result 
at z = 1 deviates about 5% from our simulation result. 
Meanwhile, ST results underpredict the simulation re- 
sults by about 10%-20% at M h > 1O 13 /i- x M . There- 
fore, the halo mass function obtained from our A^-body 
simulations supports recently proposed fitting formulae 
obtained by using large and high-resolution TV-body sim- 
ulations and vice versa. 

However, a closer look at Mh < 1O 12 /i -1 M0 in the left 
panel and Mh < 10 13 /i _1 Mq in the right panel reveals 
deviation from unity which shows that the FOF correc- 
tion may not be perfect in halos with small particles. 
Therefore, an accurate mass estimation requires keeping 
many more particles within individual halos, although 
halos with only a small number of particles (~ 20) can 
be found using the FOF algorithm. Or a more accurate 
FOF mass correction needs to be developed. In addition, 
aside from simple considerations of particle shot noise, 
there is an inherent systematic error and scatter in the 
definition of a FOF halo mass with particle number, as 
pointed out by Warren et al. [95j. 

Figure [5] is the same as Figured) but at redshift 2 and 
3 and results of the MICE fitting formula are not plotted 
because their fitting formula was fitted in the redshift 
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FIG. 4: The ratio of the halo mass function from iV-body simulations to that from three fitting formulae: Sheth and 
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panels show the results at redshift 1 and 0, respectively. The red and magenta symbols show the results with the force 
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FIG. 5: The same as Figure [4j but Crocce et al. [93] results are not plotted and the upper and bottom panels show the results 
at redshift 3 and 2. 



range 2 = 0—1. As in the results at z — and 1, our D. Halo Biases 
results are in good agreement with the Coyote results, 
although the application of their fitting formula to the 

result at z = 3 might not be reliable, because the Coyote i n Figure © we compare the halo power spectra corn- 
fitting formula was fitted between z = and z = 2. As puted from iV-body simulations with the analytical pre- 
shown in Figure H the deviation from unity at low mass dictions in real (left panel) and redshift (right panel) 
halos is also shown in Figure El space for z = 3, 2, 1, and 0.5. Halo power spectra are 

normalized by the no- wiggle fitting formula P nw (k) jll5l ] 
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FIG. 6: Comparison of halo power spectra obtained from A-body simulations to analytical predictions in real (left panel) and 
redshift (right panel) space for redshifts, z = 3, 2, 1, and 0.5. Halo power spectra are normalized by the no- wiggle fitting formula 
Pnw(fc) [1151 ] in real space and by the no- wiggle fitting formula taking the Kaiser effect into account P£ w (k) in redshift space, 
respectively. Therefore, the vertical axis shows the square of the halo bias. The solid and dotted curves show the predictions 
of the 1-loop Lagrangian resummation theory and linear theory, respectively. An enhancement of amplitude on small scales is 
due to the nonlinear effects of dynamics, bias, and redshift-space distortions. 



in real space and by the no- wiggle fitting formula taking 
the Kaiser effect into account P^ w (k) in redshift space, 
respectively. In the monopole power spectrum we have 
considered, Kaiser's enhancement effect i?hh is expressed 
as 



R 



hh 



2 1 



(31) 



from an angular average of the factor (1 + /3fi ) in 
Eq. (HH]). Therefore the vertical axis corresponds to the 
square of the halo bias, i.e., b 2 (k). The solid and dotted 
curves show the predictions of 1-loop LRT (Eqs. [S] and 
ITU) and LIN (Eqs. [TfJl and UM, respectively. When we 
calculate the prediction of the 1-loop LRT, we use the 
Coyote fit (Eq. |28|) for the halo mass function because 
the prediction of the Coyote fit is well duplicated by our 
simulation results as seen in Figures 2] and 

As increasing redshift, the amplitude of power spec- 
trum becomes larger because we impose the same halo 
mass threshold (Mh > 1.37 x 10 12 /i _1 M Q ) regardless of 
redshift. Halos which have a certain mass, for example, 
this threshold mass > 1.37x 10 12 /i _1 M , are increas- 
ingly rare as increasing redshift, and thus the amplitude 
of the power spectrum is more biased. Theoretical pre- 
dictions of 1-loop LRT and LIN well agree with A^-body 
simulations at large scales up to a certain scales. As in 
the matter power spectrum, the range of agreement in 
the 1-loop LRT prediction is wider as redshift is higher. 
Meanwhile the LIN prediction significantly deviates from 
the A^-body results as bias is higher. Therefore, the 1- 



loop LRT prediction which include higher-order bias term 
is very powerful for future high redshift BAO surveys 
to extract the cosmological information. At very large 
scales if the mass function prediction is accurate, 1-loop 
LRT and LIN predictions should reproduce the simula- 
tion result because all the higher-order Lagrangian bias 
factors (fW) should be zero except for (F'}. However, 
a closer look at small wave number k for z — 0.5 reveals 
a discrepancy in overall amplitudes between theoretical 
predictions and simulation results beyond the range of 
error bars. This might be ascribed to the fact that the 
FOF mass correction (Eq. [2^1) is not perfect in halos with 
small particles as shown in Figures HI and [S] 

In Figure [71 to examine the agreement in more quan- 
titative ways and exclude an uncertainty in amplitudes 
explained above, we eliminate dependences on the lin- 
ear bias factor and linear Kaiser's factor, dividing the 
power spectra by b 2 = (1 + (F')) 2 and b 2 R hh = (1 + 
(F')) 2 (l + 2,3/3 + /3 2 /5) in real and redshift space, re- 
spectively. For simulation results b 2 and & 2 i?hh are nu- 
merically estimated as follows. First, those factors can 
be computed among 30 realizations through 



b 2 (k) 



Phh(fc) 

Pmm {k) 



b 2 (k)R } 



hh 



-Pmm {k) 



(32) 



where P mm (k) denotes the real-space power spectrum of 
matter, and Phh(k) and P^ h {k) denote the halo power 
spectrum in real and redshift space. Those power spectra 
are obtained from A^-body simulations. Second, we use 
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FIG. 7: The same as Figured but the vertical axis is divided by b 2 = (1 + (F'}) 2 and b 2 R hh = (1 + <F')) 2 (1 + 2/3/3 + /3 2 /5) 
in real and redshift space, respectively. An enhancement of amplitude on small scales comes from the nonlinear effects of 
dynamics, bias, and redshift-space distortions. 



polynomial fitting with 

3 

b 2 {k) = A Q +Y / A t k\ (33) 
»=i 

3 

b 2 (k)R hh =B + Y,Bik\ (34) 
i=i 

Finally, computing the x 2 statistics with the data up to 
fc=0.35/iMpc _1 , we can obtain the best-fit values of scale 
independent terms, Aq and Bq, and then we substitute 
these values into b 2 and 6 2 i?hh, respectively. 

Comparing the matter power spectra in Figure [T] and 
halo power spectra in Figure the ranges of agreement 
between iV-body simulations and 1-loop LRT predictions 
with halo bias seem to be the same as those in the matter 
power spectrum for all redshift we have examined in real 
space. However, the results of the halo power spectrum 
in redshift space seem to be somewhat worse than those 
of matter power spectrum as decreasing redshift. This 
might be because the effects of nonlinear redshift-space 
distortions become larger as decreasing redshift. The 
redshift-space distortion parameter f3 ~ £l^ 55 /6 becomes 
larger as decreasing redshift since a halos at lower redshift 
are less biased than those at higher redshift. Therefore, 
1-loop LRT well replicate the iV-body simulation results 
on BAO scales even in existence of halo bias, although 
the results of low redshifts in redshift space are less accu- 
rately reproduced by 1-loop LRT. Comparing the results 
of real space at z = 0.5 and z = 1 in Figure [1] and those 
in Figure [JJ an enhancement of amplitude of halo power 
spectra on small scales is slightly moderate because (F") 



has a negative value, which is calculated from Eq. (|10l) . 
where we substitute the minimum halo mass in our sim- 
ulations, M h = 1.37 x 1O 12 /i _1 M , into Mi and infinity 
into Mi. In contrast, other redshift results [z = 2 and 
3) of halos in real space are more enhanced compared to 
those of matter power spectrum, since (F") has a posi- 
tive value at z = 2 and 3. In redshift space, amplitudes 
of the results of halo power spectrum are enhanced com- 
pared to those of matter power spectrum at all redshift. 
This is due to the nonlinear effects of bias and redshift- 
space distortions. This scale dependence of bias coming 
from the clustering of halos and nonlinear redshift-space 
distortions will be carefully examined in Figure |SJ 

Figure [5] is the same as Figure [7J but the smooth no- 
wiggle linear power spectrum P nw (fc) is replaced with 
-fmm(fc) and P^ m (k) / R mm in real and redshift space, re- 
spectively. Here P mm (k) and P^ m {k) are corresponding 
mass power spectra in real and redshift space. i? mm is 
Kaiser's enhancement factor of matter, which is obtained 
from Eq. (|3"Tj) with b = 1 as 

#mm = 1 + -/ + ~/ 2 . (35) 
3 5 

For simulation results in redshift space, b 2 Rhh/ Rmm is 
estimated from simulations through 

b 2 (k)R hh /R mm = (36) 

Using polynomial fitting with 

3 

b 2 {k)R hh /R mm = C + Y, ( 37 ) 
i=i 
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FIG. 8: The same as Figure [71 but the smooth linear power spectrum P nw (fc) is replaced with P mm (fe) and Pm m (fc)/-Rmm 
in real and redshift space, respectively in order to get rid of the nonlinearity of dynamics. Therefore deviation from unity 
shows the nonlinearity of bias and that of bias and redshift-space distortions in real and redshift space. P mm (fe) and Pm m (fc) 
are corresponding mass power spectra in real and redshift space. i? mm is Kaiser's enhancement factor of matter defined by 
R mm — 1 + 2//3 + f/5. 



we calculate the x 2 statistics with the data up to 
fc=0.35/iMpc -1 and search for the best-fit value of scale 
independent term, Cq. Then, we substitute the best- fit 
value into 6 2 i?hh/ -Rmm- 

The functional form of vertical axis in real space, 
Phh(k)/[b 2 P mm (k)], describes only the nonlinear effect of 
bias, so a deviation from unity shows the scale depen- 
dence of bias. Meanwhile in redshift space the vertical 
axis shows nonlinear effects of bias and redshift-space dis- 
tortions. We find that the scale dependence of bias does 
not show significant oscillations and are mostly smooth 
functions of scales. The scale dependence and nonlin- 
earity of bias shown in our results purely originate from 
clustering of halos themselves. LIN predicts constancy 
of halo bias on all scales. We find that scale depen- 
dences of bias calculated from our simulation results are 
in fairly good agreement with 1-loop LRT predictions up 
to k = 0.35/iMpc -1 (z = 2 and 3) within a few percent 
level and up to k = O.lh Mpc -1 (z = 2) and 0.15/iMpc -1 
(z = 3) in redshift space. Recent analytical and numer- 
ical studies claim that the scale-dependent and stochas- 
tic properties of the bias can change the redshift-space 
power spectrum and those impacts on the determin ation 
of t he growth-rate parameter would be significant (ll7l — 
Il2l| . Therefore LRT prediction is useful and powerful to 
model biasing relevant for BAO scales. 

Finally, we examine the two-point correlation function 
of halos for z — 3,2, 1, 0.5, and both in real and redshift 
space, shown in the left and right panels of Figure |H 
respectively. Two-point correlation functions of halos are 
divided by b 2 = (1 + (F'}) 2 and b 2 R hh = (1 + (F'}) 2 (1 + 



2/3/3 + /3 2 /5) in real and redshift space, respectively. The 
solid curves represent the results of 1-loop LRT while the 
dotted curves show the results of LIN. For simulation 
results b 2 and 6 2 i?hh are fitted in the same manner as in 
the case of the halo power spectrum (Figures [7] and [5J . 
Those can be calculated among 30 realizations through 

b 2 {r) = Ml, b 2 (r)R hh = J&W (38) 

smmv J Smmv' ) 

where ^ mm (r) denotes the matter correlation function 
in real space, and £hh(r) and £hh( r ) denote the halo 
correlation functions in real and redshift space. Those 
correlation functions are obtained from A-body simula- 
tions. We assume b 2 (r) and 6 2 (r)i?hh to be constant on 
sufficiently large scales both in real and redshift space 
and chose 30/i _1 Mpc < r < 80/i _1 Mpc in configuration 
space. Then we compute \ 2 f° r the above range to obtain 
b 2 and 6 2 i?hh- 

The left panel in Figure is quite similar to the left 
panel in Figure[2] This shows that halo bias does not sig- 
nificantly change the shape of BAOs peak and that the 
effects of nonlinear halo bias are not significantly relevant 
for BAO scales in the correlation function. This result 
is consistent with a recent a naly sis of halo clustering by 
using numerical simulations [116]. By comparing the left 
and right panels in Figure [5J we can see the nonlinear 
effect of redshift-space distortions. It smears the BAO 
peaks and troughs as decreasing the redshift. One-loop 
LRT with halo bias shows good agreement with A-body 
simulation results and well captures the features of non- 
linear effects as in the matter correlation function. The 
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FIG. 9: Comparison of halo two-point correlation functions computed from iV-body simulations to analytical predictions in 
real (left panel) and redshift (right panel) space for redshifts, z = 3, 2, 1, 0.5 and from bottom to top. Halo correlation 
functions are normalized by b 2 = (1 + {F')f and b 2 R hh = (1 + (F')) 2 (l + 2/3/3 + /3 2 /5) in real and redshift space, respectively. 
The solid curves represent the results of the 1-loop Lagrangian resummation theory while the dotted curves show the results 
of the linear theory. 



results of 1-loop LRT in redshift space slightly deviate 
from those of simulations at low redshifts. This should 
be related to the inaccuracy of the 1-loop LRT in the 
low-redshift power spectra as seen in Figure [7] 

V. CONCLUSION 

To exploit the full potential of upcoming high-quality 
data, we have to precisely describe observable features of 
BAOs. To achieve this purpose, taking account of the 
galaxy biasing and redshift-space distortions is essential. 
In this paper, we have used 30 large cosmological TV- 
body simulations of the standard ACDM cosmology to 
investigate the halo biases over a wide redshift range. 

First, in the matter power spectrum, 1-loop LRT is 
useful for studying nonlinear effects on scales of BAOs 
and provides better agreement with A-body results than 

1- loop SPT or the prediction of Scoccimarro [76| on large 
scales in observable redshift space. Second, in the mat- 
ter correlation function, 1-loop LRT prediction well de- 
scribes the acoustic peaks and nonlinear smearing effects 
both in real and redshift space. The 2-loop correction 
to LRT generally extends the valid range in the matter 
power spectrum. It does not have much impact on the 
correlation function in real space, because 1-loop LRT is 
already accurate enough to describe the nonlinear effects 
on BAO scales in A-body simulations [5l[. However, 
the predictions of 1-loop LRT for low-redshift correla- 
tion functions in redshift space would be improved if the 

2- loop corrections are included. 



In the halo power spectrum, we found that 1-loop 
LRT prediction has good agreement with A-body sim- 
ulation results. The ranges of agreement seem the same 
as those in the matter power spectrum for all redshift in 
real space. In redshift space, the 1-loop LRT prediction 
for the halo power spectrum is slightly worse than that 
for the matter power spectrum. This might be because 
the nonlinear effects of redshift-space distortions become 
larger as decreasing the redshift. This shows that the 
1-loop LRT prediction well reproduces the A-body sim- 
ulation results on BAO scales even in the existence of 
halo bias at least in real space. We found that the scale 
dependences of bias are pretty well reproduced by 1-loop 
LRT up to k = 0.35/iMpc- 1 (z = 2 and 3) within a few 
percent level in real space and up to k — O.l/iMpc^ 1 
(z = 2) and 0.15/iMpc^ 1 (z =3) in redshift space. In 
the halo correlation function, 1-loop LRT well describes 
nonlinear effects that smear the baryon peak and trough 
both in real and redshift space. Halo bias does not signif- 
icantly change the shape of the baryon peak. Therefore, 
the nonlinear effects of halo bias are not serious on BAO 
scales. 

Thus, LRT prediction is very powerful and reliable to 
accurately extract cosmological information for upcom- 
ing high redshift BAO surveys. 

Acknowledgments 

We greatly appreciate Takahiro Nishimichi for kindly 
providing parallelized 2nd-order Lagrangian perturba- 



14 



tion theory code. M.S. is supported by a Grant-in-Aid 
for the Japan Society for Promotion of Science (JSPS) 
fellows. T.M. acknowledges support from the Ministry 
of Education, Culture, Sports, Science, and Technol- 
ogy (MEXT), Grant-in-Aid for Scientific Research (C), 
No. 21540263, 2009. This work is supported in part by 
the JSPS Core-to-Core Program "International Research 
Network for Dark Energy" and by a Grant-in-Aid for 



Scientific Research on Priority Areas No. 467 "Probing 
the Dark Energy through an Extremely Wide and Deep 
Survey with Subaru Telescope" and by a Grant-in-Aid 
for Nagoya University Global COE Program, "Quest for 
Fundamental Principles in the Universe: from Particles 
to the Solar System and the Cosmos" , from the MEXT 
of Japan. 



[1] A. G. Ricss, A. V. Filippenko, P. Challis, A. Clocchi- 
atti, A. Diercks, P. M. Garnavich, R. L. Gilliland, C. J. 
Hogan, S. Jha, R. P. Kirshner, et al., AJ 116, 1009 
(1998). 

[2] S. Perlmutter, G. Aldering, G. Goldhaber, R. A. Knop, 
P. Nugent, P. G. Castro, S. Deustua, S. Fabbro, A. Goo- 
bar, D. E. Groom, et al., ApJ 517, 565 (1999). 

[3] C. Blake and K. Glazebrook, ApJ 594, 665 (2003). 

[4] W. Hu and Z. Haiman, Phys.Rev.D 68, 063004 (2003). 

[5] H. Seo and D. J. Eisenstein, ApJ 598, 720 (2003). 

[6] K. Glazebrook and C. Blake, ApJ 631, 1 (2005). 

[7] H. Seo and D. J. Eisenstein, ApJ 665, 14 (2007). 

[8] D. J. Eisenstein, W. Hu, and M. Tegmark, ApJ 504, 
L57+ (1998). 

[9] A. Cooray, W. Hu, D. Huterer, and M. Joffre, ApJ 557, 
L7 (2001). 

[10] W. Hu, in Observing Dark Energy, edited by S. C. Wolff 
& T. R. Lauer (2005), vol. 339 of Astronomical Society 
of the Pacific Conference Series, pp. 215 — h 

[11] E. V. Linder, Phys.Rev.D 68, 083504 (2003). 

[12] T. Matsubara, ApJ 615, 573 (2004). 

[13] L. Amendola, C. Quercellini, and E. Giallongo, MNRAS 
357, 429 (2005). 

[14] A. Albrecht, G. Bernstein, R. Cahn, W. L. Freedman, 
J. Hewitt, W. Hu, J. Huth, M. Kamionkowski, E. W. 
Kolb, L. Knox, et al., arXiv:astro-ph/0609591 (2006). 

[15] D. Dolney, B. Jain, and M. Takada, MNRAS 366, 884 
(2006). 

[16] A. Albrecht, L. Amendola, G. Bernstein, D. Clowe, 
D. Eisenstein, L. Guzzo, C. Hirata, D. Huterer, R. Kir- 
shner, E. Kolb, et al., arXiv:0901.0721 (2009). 

[17] D. J. Schlegel, C. Bebek, H. Heetderks, S. Ho, M. Lamp- 
ton, M. Levi, N. Mostek, N. Padmanabhan, S. Perlmut- 
ter, N. Roe, et al., arXiv:0904.0468 (2009). 

[18] A. Refregier, A. Amara, T. D. Kitching, A. Rassat, 
R. Scaramella, J. Weller, and f. t. Euclid Imaging Con- 
sortium, arXiv:1001.0061 (2010). 

[19] S. Cole, W. J. Percival, J. A. Peacock, P. Norberg, C. M. 
Baugh, C. S. Frenk, I. Baldry, J. Bland-Hawthorn, 
T. Bridges, R. Cannon, et al., MNRAS 362, 505 (2005). 

[20] D. J. Eisenstein, I. Zehavi, D. W. Hogg, R. Scoccimarro, 
M. R. Blanton, R. C. Nichol, R. Scranton, H. Seo, 
M. Tegmark, Z. Zheng, et al., ApJ 633, 560 (2005). 

[21] M. Tegmark, D. J. Eisenstein, M. A. Strauss, D. H. 
Weinberg, M. R. Blanton, J. A. Frieman, M. Fukugita, 
J. E. Gunn, A. J. S. Hamilton, G. R. Knapp, et al., 
Phys.Rev.D 74, 123507 (2006). 

[22] N. Padmanabhan, D. J. Schlegel, U. Seljak, A. Makarov, 
N. A. Bahcall, M. R. Blanton, J. Brinkmann, D. J. 
Eisenstein, D. P. Finkbeiner, J. E. Gunn, et al., MN- 
RAS 378, 852 (2007). 



[23] W. J. Percival, R. C. Nichol, D. J. Eisenstein, J. A. 
Frieman, M. Fukugita, J. Loveday, A. C. Pope, D. P. 
Schneider, A. S. Szalay, M. Tegmark, et al., ApJ 657, 
645 (2007). 

[24] T. Okumura, T. Matsubara, D. J. Eisenstein, I. Kayo, 

C. Hikage, A. S. Szalay, and D. P. Schneider, ApJ 676, 
889 (2008). 

[25] B. A. Reid, W. J. Percival, D. J. Eisenstein, L. Verde, 

D. N. Spergel, R. A. Skibba, N. A. Bahcall, T. Budavari, 
J. A. Frieman, M. Fukugita, et al., MNRAS 404, 60 
(2010). 

[26] N. Kaiser, ApJ 284, L9 (1984). 

[27] J. M. Bardeen, J. R. Bond, N. Kaiser, and A. S. Szalay, 

ApJ 304, 15 (1986). 
[28] N. Kaiser, MNRAS 227, 1 (1987). 

[29] A. Meiksin, M. White, and J. A. Peacock, MNRAS 304, 
851 (1999). 

[30] H.-J. Seo and D. J. Eisenstein, ApJ 633, 575 (2005). 
[31] D. J. Eisenstein, H.-J. Seo, and M. White, ApJ 664, 
660 (2007). 

[32] E. Huff, A. E. Schulz, M. White, D. J. Schlegel, and 

M. S. Warren, Astroparticle Physics 26, 351 (2007). 
[33] H.-J. Seo, E. R. Siegel, D. J. Eisenstein, and M. White, 

ApJ 686, 13 (2008). 
[34] R. E. Smith, R. Scoccimarro, and R. K. Sheth, 

Phys.Rev.D 77, 043525 (2008). 
[35] N. Padmanabhan and M. White, Phys.Rev.D 80, 

063508 (2009). 

[36] H.-J. Seo, J. Eckel, D. J. Eisenstein, K. Mehta, 
M. Metchnik, N. Padmanabhan, P. Pinto, R. Takahashi, 
M. White, and X. Xu, ApJ 720, 1650 (2010). 

[37] K. T. Mehta, H.-J. Seo, J. Eckel, D. J. Eisenstein, 
M. Metchnik, P. Pinto, and X. Xu, ApJ 734, 94 (2011). 

[38] M. R. Blanton, D. Eisenstein, D. W. Hogg, and I. Ze- 
havi, ApJ 645, 977 (2006). 

[39] R. E. Smith, R. Scoccimarro, and R. K. Sheth, 
Phys.Rev.D 75, 063512 (2007). 

[40] A. G. Sanchez and S. Cole, MNRAS 385, 830 (2008). 

[41] C. Blake, S. Brough, M. Colless, C. Contreras, 
W. Couch, S. Croom, T. Davis, M. J. Drinkwater, 
K. Forster, D. Gilbank, et al., arXiv: 1104.2948 (2011). 

[42] H. J. Mo and S. D. M. White, MNRAS 282, 347 (1996). 

[43] J. A. Peacock and R. E. Smith, MNRAS 318, 1144 
(2000). 

[44] U. Seljak, MNRAS 318, 203 (2000). 

[45] A. Cooray and R. Sheth, Phys. Rep. 372, 1 (2002). 

[46] R. E. Smith, R. K. Sheth, and R. Scoccimarro, 

Phys.Rev.D 78, 023523 (2008). 
[47] A. Elia, S. Kulkarni, C. Porciani, M. Pietroni, and 

S. Matarrese, arXiv:1012.4833 (2010). 
[48] T. Matsubara, Phys.Rev.D 83, 083518 (2011). 



15 



[49] T. Matsubara, Phys.Rev.D 78, 083519 (2008). 

[50] F. Bernardeau, S. Colombi, E. Gaztanaga, and R. Scoc- 

cimarro, Phys. Rep. 367, 1 (2002). 
[51] T. Okamura, A. Taruya, and T. Matsubara, 

arXiv: 1105. 1491 (2011). 
[52] T. Padmanabhan, Structure Formation in the Universe 

(1993). 

[53] T. Buchert, A&A 223, 9 (1989). 

[54] T. Buchert, MNRAS 254, 729 (1992). 

[55] T. Buchert and J. Ehlers, MNRAS 264, 375 (1993). 

[56] T. Buchert, MNRAS 267, 811 (1994). 

[57] P. Catelan, MNRAS 276, 115 (1995). 

[58] E. Hivon, F. R. Bouchet, S. Colombi, and 
R. Juszkiewicz, A&A 298, 643 (1995). 

[59] J. Ehlers and T. Buchert, General Relativity and Grav- 
itation 29, 733 (1997). 

[60] T. Tatekawa, arXiv:astro-ph/0412025 (2004). 

[61] T. Matsubara, Phys.Rev.D 77, 063530 (2008). 

[62] Y. Suto and M. Sasaki, Physical Review Letters 66, 264 

(1991) . 

[63] N. Makino, M. Sasaki, and Y. Suto, Phys.Rev.D 46, 585 

(1992) . 

[64] B. Jain and E. Bertschinger, ApJ 431, 495 (1994). 
[65] T. Nishimichi, H. Ohmuro, M. Nakamichi, A. Taruya, 

K. Yahata, A. Shirata, S. Saito, H. Nomura, K. Ya- 

mamoto, and Y. Suto, PAS J 59, 1049 (2007). 
[66] J. Carlson, M. White, and N. Padmanabhan, 

Phys.Rev.D 80, 043531 (2009). 
[67] M. Crocce and R. Scoccimarro, Phys.Rev.D 73, 063519 

(2006). 

[68] M. Crocce and R. Scoccimarro, Phys.Rev.D 73, 063520 
(2006). 

[69] M. Crocce and R. Scoccimarro, Phys.Rev.D 77, 023533 
(2008). 

[70] P. Valageas, A&A 465, 725 (2007). 
[71] M. Pietroni, JCAP 10, 36 (2008). 

[72] S. Matarrese and M. Pietroni, Modern Physics Letters 

A 23, 25 (2008). 
[73] A. Taruya and T. Hiramatsu, ApJ 674, 617 (2008). 
[74] J. C. Jackson, MNRAS 156, IP (1972). 
[75] W. L. W. Sargent and E. L. Turner, ApJ 212, L3 (1977). 
[76] R. Scoccimarro, Phys.Rev.D 70, 083007 (2004). 
[77] A. J. S. Hamilton, in The Evolving Universe, edited by 

D. Hamilton (1998), vol. 231 of Astrophysics and Space 

Science Library, pp. 185 — K 
[78] J. L. Tinker, D. H. Weinberg, and Z. Zheng, MNRAS 

368, 85 (2006). 
[79] A. Taruya, T. Nishimichi, S. Saito, and T. Hiramatsu, 

Phys.Rev.D 80, 123503 (2009). 
[80] A. Taruya, T. Nishimichi, and S. Saito, Phys.Rev.D 82, 

063522 (2010). 

[81] E. A. Kazin, A. G. Sanchez, and M. R. Blanton, 

arXiv:1105.2037 (2011). 
[82] J. Kwan, G. F. Lewis, and E. V. Linder, 

arXiv:1105.1194 (2011). 
[83] B. A. Reid and M. White, arXiv: 1105.4165 (2011). 
[84] J. Tang, I. Kayo, and M. Takada, arXiv:1103.3614 

(2011). 

[85] A. Taruya, S. Saito, and T. Nishimichi, Phys.Rev.D 83, 

103527 (2011). 
[86] C. Lacey and S. Cole, MNRAS 262, 627 (1993). 
[87] V. R. Eke, S. Cole, and C. S. Frenk, MNRAS 282, 263 

(1996). 

[88] T. T. Nakamura and Y. Suto, Progress of Theoretical 



Physics 97, 49 (1997). 
[89] J. P. Henry, ApJ 534, 565 (2000). 

[90] A. Jenkins, C. S. Frenk, S. D. M. White, J. M. Col- 
berg, S. Cole, A. E. Evrard, H. M. P. Couchman, and 
N. Yoshida, MNRAS 321, 372 (2001). 

[91] W. H. Press and P. Schechter, ApJ 187, 425 (1974). 

[92] J. R. Bond, S. Cole, G. Efstathiou, and N. Kaiser, ApJ 
379, 440 (1991). 

[93] R. K. Sheth, H. J. Mo, and G. Tormen, MNRAS 323, 
1 (2001). 

[94] R. K. Sheth and G. Tormen, MNRAS 308, 119 (1999). 
[95] M. S. Warren, K. Abazajian, D. E. Holz, and 

L. Teodoro, ApJ 646, 881 (2006). 
[96] D. S. Reed, R. Bower, C. S. Frenk, A. Jenkins, and 

T. Theuns, MNRAS 374, 2 (2007). 
[97] M. Crocce, P. Fosalba, F. J. Castander, and 

E. Gaztanaga, MNRAS 403, 1353 (2010). 
[98] M. Manera, R. K. Sheth, and R. Scoccimarro, MNRAS 

402, 589 (2010). 
[99] S. Bhattacharya, K. Heitmann, M. White, Z. Lukic, 

C. Wagner, and S. Habib, ApJ 732, 122 (2011). 
[100] H. J. Mo, Y. P. Jing, and S. D. M. White, MNRAS 284, 

189 (1997). 

[101] V. Springel, MNRAS 364, 1105 (2005). 

[102] E. Komatsu, K. M. Smith, J. Dunkley, C. L. Bennett, 
B. Gold, G. Hinshaw, N. Jarosik, D. Larson, M. R. 
Nolta, L. Page, et al., ApJS 192, 18 (2011). 

[103] M. Crocce, S. Pueblas, and R. Scoccimarro, MNRAS 
373, 369 (2006). 

[104] T. Nishimichi, A. Shirata, A. Taruya, K. Yahata, 
S. Saito, Y. Suto, R. Takahashi, N. Yoshida, T. Mat- 
subara, N. Sugiyama, et al., PASJ 61, 321 (2009). 

[105] A. Lewis, A. Challinor, and A. Lasenby, ApJ 538, 473 
(2000). 

[106] P. Valageas and T. Nishimichi, A&A 527, A87+ (2011). 

[107] M. Davis, G. Efstathiou, C. S. Frenk, and S. D. M. 
White, ApJ 292, 371 (1985). 

[108] R. W. Hockney and J. W. Eastwood, Computer Sim- 
ulation Using Particles (New York: Taylor & Francis 
Group, 1988). 

[109] R. E. Angulo, C. M. Baugh, C. S. Frenk, and C. G. 
Lacey, MNRAS 383, 755 (2008). 

[110] R. Takahashi, N. Yoshida, M. Takada, T. Matsubara, 
N. Sugiyama, I. Kayo, A. J. Nishizawa, T. Nishimichi, 
S. Saito, and A. Taruya, ApJ 700, 479 (2009). 

[Ill] Y. P. Jing, ApJ 620, 559 (2005). 

[112] R. Takahashi, N. Yoshida, T. Matsubara, N. Sugiyama, 
I. Kayo, T. Nishimichi, A. Shirata, A. Taruya, S. Saito, 
K. Yahata, et al., MNRAS 389, 1675 (2008). 

[113] U. Seljak, N. Hamaus, and V. Desjacques, Physical Re- 
view Letters 103, 091303 (2009). 

[114] N. Hamaus, U. Seljak, V. Desjacques, R. E. Smith, and 
T. Baldauf, Phys.Rev.D 82, 043515 (2010). 

[115] D. J. Eisenstein and W. Hu, ApJ 496, 605 (1998). 

[116] A. G. Sanchez, C. M. Baugh, and R. Angulo, MNRAS 
390, 1470 (2008). 

[117] D. Jeong and E. Komatsu, ApJ 691, 569 (2009). 

[118] V. Desjacques and R. K. Sheth, Phys.Rev.D 81, 023526 
(2010). 

[119] E. Jennings, C. M. Baugh, and S. Pascoli, MNRAS 410, 
2081 (2011). 

[120] T. Okumura and Y. P. Jing, ApJ 726, 5 (2011). 
[121] S. Saito, M. Takada, and A. Taruya, Phys.Rev.D 83, 
043529 (2011). 



